A PRELIMINARY STUDY FOR A 
TETRAHEDRON FORMATION: 
QUALITY FACTORS AND VISUALIZATION 

Jose J. Guzman* 

Conrad Schif? 

a A. solutions, Inc. 

Mission Design Division 
10001 Derekwood Lane , Suite 215 
Lanham , MD 20706 , USA 

Abstract 

Spacecraft flying in tetrahedron formations are excellent for electromagnetic and plasma studies. 

The quality of the science recorded is strongly affected by the tetrahedron evolution. This paper ls 
a preliminary study on the computation of quality factors and visualization for a formation of four 
or five satellites. Four of the satellites are arranged geometrically in a tetrahedron shape. If a hith 
satellite is present, it is arbitrarily initialized at the geometric center of the tetrahedron. The fifth 
satellite could act as a collector or as a spare spacecraft. Tetrahedron natural coordinates are employed 
for the initialization. The natural orbit evolution is visualized in geocentric equatorial inertial and in 
geocentric solar magnetospheric coordinates. 


INTRODUCTION 

Spacecraft formation flying and/or distributed space 
systems allow measurements not possible with single 
spacecraft. For space science missions, observable field 
parameters (particle populations, electric and mag- 
netic fields) vary both in space and time. 1 Thus, 
understanding of the processes within the field re- 
quires accurate and precise observations and analysis 
of both the temporal and spatial variations. A tetra- 
hedron formation is utilized since four spacecraft, with 
adjustable separations, are the minimum needed to 
resolve a three-dimensional structure, at least to the 
lowest order in the physical field gradients. 2 

The appropriate spacecraft separations, of course, 
depend on the field or event being sampled. As a 
result, formation flying missions require closer inter- 
action between the project scientists and mission de- 
signers than a typical single-spacecraft mission. A 
means of communicating quantitatively is to define 
quality factors that affect the data collection and anal- 
ysis processes. These quality factors are shape- and 
size-dependent parameters that are functions of the 
geometric shape and evolution of the formation. Once 
the scientists have defined the quality factors appro- 
priate for their mission, the mission analysts can try 
to tailor the trajectory design appropriately. Further- 
more, to communicate qualitatively, visualization of 
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the formation — in coordinate frames useful to the 
scientists is important. This paper is a preliminary 
study on the computation of quality factors and visu- 
alization for a formation of four or five satellites. 

TETRAHEDRON MISSIONS 

Several missions are planning to use tetrahedron for- 
mations for field and plasma studies. In fact, one such 
mission, Cluster IT, is currently flying and operat- 
ing successfully. 3 The Cluster II mission was launched 
in pairs using Russian Soyuz rockets: two spacecraft 
launched on July 12, 2000 and two on August 9, 
2000. Cluster's primary mission is to explore the Sun- 
Earth connection, specifically the interaction of the 
Sun and the Earth’s magnetosphere. In the process 
of planning the Cluster mission, fuel optimization of 
the maneuvers needed to initialize, modify and main- 
tain the formation was considered by J. Rodriguez- 
Canabal and M. Bello-Mora. 4 Later a different opti- 
mization method was employed by J. Schoenmaekers. 5 
In J. Schoenmaekcrs’s paper, both methods are com- 
pared for one mission scenario. Moreover, a further- 
developed strategy is presented by M. Bello-Mora and 
Rodriguez-Canabal. 6 Implementation and operational 
results for the maneuvers are reported by D. Hockens 
and J. Schoenmaekers. 7 

* Cluster II is a replacement of the original Cluster mission 
which was lost in a launch failure during the maiden flight of 
the Ariane 5 rocket on June 4, 1996. 
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The Auroral Lites mission concept also proposed us- 
ing four spacecraft. This mission’s main goal would 
have been to acquire an unprecedented understanding 
of the multiscale dynamics of space plasmas by ex- 
ploring the Earth’s auroral zone . 8,9 Auroral Lites also 
would have been a technology pathfinder for flying for- 
mations of small spacecraft. Unfortunately, the mis- 
sion did not make the final selection of the 1998 NASA 
Medium-class Explorers (MIDEX) program. An- 
other proposed mission that will utilize a tetrahedron 
formation is the Magnetospheric Multiscale (MMS) 
Mission . 10 MMS will determine the small-scale basic 
plasma processes which transport, accelerate and en- 
ergize plasmas in thin boundary and current layers. 
These processes control the structure and dynamics of 
the Earth’s magnetosphere. Preliminary mission de- 
sign and analysis for the MMS reference path, which 
includes a double lunar swingby for one of its phases, 
has been presented by A. Edery and C. Schiff. 

THE VOLUMETRIC TENSOR 

The spatial gradient of a field observed with multiple 
spacecraft is dependent on the inverse of a symmet- 
ric tensor formed from the positions of the spacecraft 
(see Harvey 1 ). Harvey adds that “the importance of 
the volumetric tensor for describing the geometry of 
a tetrahedron was first noted by J. Schoenmaekers 
of the European Space Operations Centre (ESOC) 
Flight Dynamics Division”. Thus, this tensor is impor- 
tant for both the gradient calculations and for providing 
information about the shape and orientation of the for- 
mation . For N spacecraft, in an arbitrary frame I and 
about the mesocenter (mean position or centroid), 
the volumetric tensor can be expressed as, 

(1) 

i=i 

where the double bar is used to indicate a second order 
tensor, and r mi is the position (column) vector of the i~ 
th spacecraft relative to the mesocenter, me. Relative 
to the origin of the arbitrary frame /, the mesocenter 
is computed as, 

= JfY.fi, (2) 

i— 1 

where r* is the position vector of the i-th spacecraft 
relative to the origin of the arbitrary frame /. The 
mesocenter is equal to the system’s center of mass 
if all the spacecraft have the same mass. Furthermore, 
if the tetrahedron is regular (the separations between 
vertices are equal), the mesocenter equals the geo- 
metric center (a point equidistant from all vertices). 
For a tetrahedron, it can be shown that the volume V 
is given by, 

V = »^\, ( 3 ) 

3 


where \R I/mc \ is the determinant of the tensor. Hence 
the name, volumetric tensor. 

R elationship to the Inertia Tensor 

The inertia tensor about the mesocenter is closely 
related to the volumetric tensor. In fact, the eigenvec- 
tors are the same and the eigenvalue associated with 
the largest (smallest) eigenvalue of the inertia tensor 
is associated with the smallest (largest) eigenvalue of 
the volumetric tensor. See Harvey 1 for the equation 
relating the eigenvalues. 

Pseudo-Ellipsoid 

Analogous to the inertia-ellipsoid, a pseudo-ellipsoid 
can be constructed with the eigenvalues and eigenvec- 
tors of the volumetric tensor. Therefore, the pseudo- 
ellipsoid has information about the shape and orienta- 
tion of the formation. This proves very useful in the 
formulation of quality factors. 


NATURAL COORDINATES 

For both mission design and science purposes, it 
might be convenient to sometimes work in a coordinate 
frame that is dependent on the shape of the tetrahe- 
dron. This can be accomplished by utilizing natural 
coordinates. Natural coordinates are typically used in 
finite element analysis . 12 Analogous to the volumetric 
tensor, the tetrahedron natural coordinates are depen- 
dent on the volume of the tetrahedron. A tetrahedron 
in rectangular and natural coordinates is displayed in 

Figure 1. The faces of the tetrahedron are numbered 

S /C 4 

Face 2 (fr = 0 ) 


Face 3 (£ 



Fig. 1 Tetrahedron Natural Coordinates 

such that the i-th face is opposite to the z-th vertex 
(spacecraft). The tetrahedral coordinates for an in- 
terior point, £ = {£i, £ 2 j £ 3 * £ 4 } T > are defined as the 
ratio of the volume of the subtetrahedron formed by 
the interior point and face i to the tetrahedron volume. 
That is, £ = {Vi/V, V 2 /V, \ WE, V 4 /V} t . Thus, the 
i-th face has £* = 0. Furthermore, note that 


Vi + V2 + V3 + V4 — V , ( 4 ) 


which implies that, 

£1 + £2 + £3 + £4 — L 

Therefore, with this last equation and relating the 
cartesian coordinates {x, y, 2 } of an interior point to 
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Fig. 2 Example: Inertial Coordinates 


the tetrahedral coordinates, 
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( 6 ) 


where the cartesian coordinates of the points that de- 
fine the tetrahedron are given by the vectors f x = 
{x*, yi , Zi} T for i- 1, ... ,4. For example, in natural 
coordinates the geometric center of a regular tetrahe- 
dron is located at £ — {1/4, 1/4, 1/4, 1/4} 7 . Then, 
using Equation (6) the geometric center cartesian co- 
ordinates are readily available. See Figure 2 for an 
example with orbital data. The transformation from 
cartesian to natural coordinates is given by taking the 
inverse of the matrix in Equation (6), provided the 
tetrahedron volume is non-zero. 


in defining the formation geometry and stationkeep- 
ing requirements. Before introducing some of them, 
various parameter definitions are required. 

Some Geometric Definitions 

First, recall that in a regular tetrahedron the sep- 
arations between vertices are equal. Then, the ideal 
volume and surface are defined as the volume and 
surface area of a regular tetrahedron with vertex sepa- 
rations equal to the average of the 6 distances between 
the true, or actual, vertices. The circumscribing 
sphere is defined as the sphere with center at the ge- 
ometric center and radius equal to the distance from 
that center to any spacecraft (thus, ail 4 spacecraft lie 
on its surface). 

One- Dimensional Parameters 

The following parameters are useful for comparing 
the true tetrahedron against an ideal (regular) tetra- 
hedron. A regular tetrahedron might be useful in cases 
where sampling of the structure of a field is not as im- 
portant as understanding its transient or fluctuating 
events. 13 

• The Glassmeier Q G M parameter is defined as, 

True Volume True Surface . . 

O — -I — 4- 1. (7) 

Ideal Volume Ideal Surface 

This parameter takes values between 1 (the four 
spacecraft are in a line) and 3 (a regular tetrahe- 
dron). 

• The Robert/Roux Q RR parameter is defined as, 

( 9k True Volume \ ^ 

Qrr - ^ 2^/3 Sphere Volume / 


Aside: Volume Integrals 

It should be noted that natural coordinates simplify 
the integration over the volume of a tetrahedron. This 
may be useful for some of the science analysis. In fact, 
the integration limits for each variable become con- 
stants. For example, for an integral, k, of a function, 
/, over the volume of the tetrahedron, 

k = f f(x,y,z)dV 
J Vo 1 

r i ri-u r 

= / / / /«i, 6. €».€») <*&<*&<*&• 

Jo Jo Jo 

As a result, the integrations may be easily obtained 
numerically and even exactly. See Yang 12 for details 
and examples. 

QUALITY FACTORS 

The quality factors are shape- and size-dependent 
parameters that affect the resolution of spatial gradi- 
ents. Many have been examined and compared, but 
in this paper the ones examined by P. Robert et al. 13 
are utilized. In general, the quality factors can aid 


where the sphere volume is the volume of the 
circumscribing sphere. This parameter has a min- 
imum value of 0 and equals 1 for a regular tetra- 
hedron. 


• The Q sr parameter is defined as, 


Qsr 


1 

2 


( a+6+c 
5 



( 9 ) 


This parameter uses information from the volu- 
metric tensor, namely its eigenvalues. Specifically, 
referring to the volumetric tensor, let a, b , c be its 
largest, middle and smallest eigenvalues respec- 
tively. 

• The Q R8 parameter is defined as, 

_ True Volume qqn 

^ R8 Ideal Volume* 


This parameter is the first term of Equation (7). It 
has a minimum value of zero and equals 1 for a regular 
tetrahedron. 
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Two-Dimensional Parameter Space 


APPROACH 


Now, in an effort to have some parameters that 
describe the shape of the formation better, a two- 
dimensional parameter space can be defined by utiliz- 
ing the concepts of elongation (prolateness) and pla- 
narity (oblateness) of the pseudo-ellipsoid. Following 
P. Robert et al., 13 the elongation and planarity are 
computed as, 


E = 1 -(6/a), (11) 

P = 1 — (c/6). (12) 

Different characteristic shapes can de defined in terms 
of the values of E and P. Both parameters take val- 
ues between 0 and 1. In fact, if E = P — 0, the 
tetrahedron is regular; if E = 1 the four spacecraft 
are co-linear; and, if P = 1 the four spacecraft are 
co-planar. 13 


SCIENTIFIC COORDINATE FRAMES 

Several coordinate frames are useful for electro- 
magnetic and plasma studies. 14 Thus, these frames 
are important for visualization of the satellite forma- 
tion in a frame relevant to the scientists. For now, two 
frames are considered: (1) the Solar Magnetospheric 
(GSM) coordinate system, and (2) the Solar Magnetic 
(SM) coordinate system. These frames are defined as 
follows. First, let r 0G he a unit vector in the Earth to 
Sun direction* and rid he a unit vector in the magnetic 
north direction.* Then, the two frames of interest are 
defined, in terms of the axis directions, as 

L GSM coordinate system: X GSM = f e Q , T G sm = 
Ilflj Xfgsoll ’ ZGSM _ AgSM r ° SM ’ ^ 


2. SM coordinate system: Z SM = hd, Ui 

ndxr^o v V" ^7 

Asm ” Ysm x Zsm ' 


In both cases, the Earth dipole axis is in the X-Z 
plane. The angle between the north magnetic pole, 
fid , and the Z G SM axis is called the dipole tilt angle. 
Thus, the difference between these two frames is a ro- 
tation about the Y axes (which are the same for both 
frames). The rotation angle is simply the dipole tilt 
angle. In this paper, only the GSM tetrahedron evo- 
lution is shown. 


f The Earth to Sun vector is computed using either ephemeris 
or an approximation. 

*The dipole axis is computed using dipole angles from the 
International Geomagnetic Reference Field (IGRF). The IGRF 
is a series of mathematical models of the Earth’s main magnetic 
field and its secular variation. Each model comprises a set of 
spherical harmonic (or Gauss) coefficients in a series expansion 
of the geomagnetic potential . 15 


To test all the concepts introduced, several runs 
were performed. Specifically, the quality factors are 
computed for the tetrahedrons in Table 1. The pa- 
rameters associated with these test cases are discussed 
next. 


Formation Number 

Tetrahedron Spacecraft 

1 

(1.2,3, 4) 

2 

(1,2, 3, 5) 

3 

(1,2, 4, 5) 

4 

(1.3, 4,5) 

5 

(2, 3, 4, 5) 


Table 1 Tetrahedron Permutations 


Software 

For this study several pieces of software are uti- 
lized: FreeFlyer®, 16 MATLAB® and some internal 
legacy C code that transforms to and from the inertial 
and scientific frames. FreeFlyer ® has an embedded 
scripting language with many capabilities, including 
calling MATLAB ® and other programs (calls can be 
made in both directions). In this case, the legacy 
C code is compiled as a MATLAB® MEX file that 
FreeFlyer® can then access. 

Force Model 

All spacecraft are assumed to have the same physi- 
cal properties: a dry mass of 200 kg, an area of 0.6613 
m 2 , a drag coefficient of 2.2, and a 1.4 reflectivity 
coefficient. The FreeFlyer® integrator selected uses 
the Runge-Kutta-Verner 8(9) formulas. 17 The force 
model selected has the following contributions: 8x8 
Earth potential model, atmospheric drag, solar radi- 
ation pressure, and the Sun and the Moon as point 
masses. 

Reference Orbit 

The reference orbit is selected to be the same as in 
Hametz et al. 9 For scientific reasons, it is designed as a 
500 x 7000 km orbit at an Earth equatorial inclination 
of 80 degrees. The elements used to define the initial 
state are in J 2 Brouwer-Lyddane elements (BLJ2). 

Initializing the Formation 

A regular tetrahedron is utilized as the shape re- 
quired to meet the science objectives. An initial 
separation of 10 km is used. Recall that a regu- 
lar tetrahedron might be useful in cases where un- 
derstanding transient and fluctuating events 13 is im- 
portant. The initial formation geometry is specified 
via the FreeFlyer® close formation utility. In this 
utility, a reference orbit is specified and the individ- 
ual spacecraft locations are specified in the (rotating) 
velocity-binormal-normal (VBN) frame. Specifically, 
each spacecraft location is specified using spherical co- 
ordinates in the VBN frame. 18 In terms of velocity, 
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the velocity direction is kept the same as that of the 
reference orbit, i.e. the VBN frame velocity direc- 
tion. Then, the velocity magnitude is obtained using 
the tw o-body proble m (2BP) energy equation, thus, 
v = y/fi(2/r - 1/a), where v is the velocity magni- 
tude, r is the position magnitude (radial distance), p 
is the gravitational parameter and a is the semi-major 
axis. 

The fifth spacecraft is initialized to be at the geo- 
metric center of the formation. (This is easily done 
by utilizing tetrahedron natural coordinates). It is as- 
sumed that these initial conditions can be achieved by 
the spacecraft from some initial orbit (s). Of course , 
different initial conditions should be tried to lower the 
fuel cost to initialize the formation . This optimization 
process will depend on the pre-initialization orbital 
states which in turn depend on the launch to orbit 
strategy. Examples of formation initialization are dis- 
cussed by Mailhe et al. 19 and by Hughes and Mailhe. 20 

A Note on Stationkeeping and Relative Orbits 

In this paper the formation stationkeeping and/or 
re-configuration problem is not considered. Never- 
theless, many efforts are in progress in this arena for 
different types of missions. See for instance the work 
by Chao et al. 21 In terms of relative orbit design, the 
paper by Alfriend and Schaub provides a good sum- 
mary and some possible solutions for Earth orbiting 
platforms. 22 It should be noted that for close for- 
mations and/or precise formations, relative navigation 
and communication is most likely needed. 

RESULTS 

Plots of quality factors as a function of BLJ2 mean 
anomaly for one orbit and for the 5 tetrahedron per- 
mutations (Table 1) are displayed in Figures 3-6. The 
tetrahedron is initialized at 4 different points corre- 
sponding to mean anomalies of 0 (Figure 3), 90 (Figure 
4), 180 (Figure 5) and 270 degrees (Figure 6). It ap- 
pears that switching from one formation permutation 
to another might allow the maximization of a particu- 
lar quality factor. This would allow taking advantage 
of the natural dynamics and should be further in- 
vestigated. In general, the location along the orbit 
chosen for the formation initialization should represent 
the region where the primary science will be carried 
out. Thus, for space science missions, it is possible 
to have two formation control regimes. The first one 
would be a “no-science control” regime where the goal 
is to maintain the spacecraft within a certain maxi- 
mum separation that will facilitate a smooth transition 
into the “science control” regime where the formation 
geometry and/or quality factors is enforced to the re- 
quired (and feasible) tolerances. 

Keplerian Turning Points 

In most of the quality factor figures, note the pres- 
ence of at least 5 turning points. These points take 


place at true anomalies of 0, 90, 180, 270 and back to 
0. (The plots are with respect to BLJ2 Mean Anomaly, 
thus, in this case, the points are at mean anomalies of 
approximately 0, 54, 180, 306 and back to 0). The 
presence of the turning points is not surprising since, 
from the 2BP, a Keplerian orbit has a radial distance 
rate given by 

, _ eh sintfl*) , (13) 

P 

where e is the orbit eccentricity, h = f x r is the spe- 
cific angular momentum, 0* is the true anomaly and 
p = a(l — e 2 ) is the semi-latus rectum. Thus, the max- 
imum and minimum radial rate changes occur at true 
anomalies of 90 and 270 degrees, respectively; and the 
radial rate is 0 at the apses (true anomalies of 0 and 
180 degrees). 

The geometry of the formation (initialized at a BLJ2 
mean anomaly of 180 degrees or at apogee) is shown 
at the 5 turning points in Figures 7 and 8 in the GSM 
frame. For this particular initialization, note that for 
true anomalies of 270 and 90 degrees, the 4th and 5th 
spacecraft appear to be on top of each other! In fact, 
the distance between these spacecraft at these two 
points is about 14 and 20 meters respectively. This 
case, therefore, would need some type of active con- 
trol to maintain the 5th spacecraft at a safe distance. 

SUMMARY 

This study was an effort to compute the quality 
factors associated with tetrahedron formations and 
to visualize the geometric evolution of the formation. 
Tetrahedron permutations were examined by consid- 
ering a formation of four satellites and adding a fifth 
satellite. The fifth satellite was placed at the geometric 
center of the tetrahedron chosen for the first four satel- 
lites. The concept of the volumetric tensor and quality 
factors were reviewed and implemented. Also, tetrahe- 
dron natural coordinates were introduced. These coor- 
dinates might be very useful for further studies. Other 
initial conditions for the fifth satellite and switching 
from one formation permutation to another in order to 
maximize a particular quality factor should be investi- 
gated. In general, the location along the orbit chosen 
for the formation initialization should represent the 
region where the primary science will be carried out. 
Of course, this location will affect the quality factors. 
Also, varying the initial distances in the formation 
should be studied. It is of interest to know how the sci- 
ence is affected by the evolution of the tetrahedron, a 
key factor in defining the stationkeeping requirements. 
Furthermore, the quality factors could be used in a 
feedback mechanism to control the formation shape. 

ACKNOWLEDGEMENTS 

Part of this work was done under NASA contract 
NASS-01090/Task 115. Special thanks to Dr. Steven 


American Institute of Aeronautics and Astronautics 


A. Curtis and Mr. David C. Folta from NASA God- 
dard Space Flight Center. 

REFERENCES 

l C. Harvey, “Spatial Gradients and the Volumetric Tensor,” 
Analysis Methods for Multi- Spacecraft Data (G. Paschmann and 
P. Daly, eds.), pp. 307-322, Noordwijk, The Netherlands: ISSI 
Report SR-001, ESA Publications Division, 1998. 

2 G. Paschmann and P. Daly, “Introduction,” Analysis 
Methods for Multi- Spacecraft Data (G. Paschmann and P. Daly, 
eds.), Noordwijk, The Netherlands: ISSI Report SR-001, ESA 
Publications Division, 1998. 

3 R. Mugellesi-Dow, J. Dow, and G. Gienger, “Cluster II: 
Prom Launch to First Constellation,” 16th International Sym- 
posium on Space Flight Dynamics , Pasadena, California, De- 
cember 2001. 

4 J. Rodriguez-Canabal and M. Bello-Mora, “Cluster: Con- 
solidated Report on Mission Analysis,” European Space Opera- 
tion Centre, Technical Report, Darmstadt, Germany, July 1990. 
ESOC CL-ESC-RP-0001. 

5 J. Schoenmaekers, “Cluster: Fuel Optimum Spacecraft 

Formation Control,” ESA Symposium on Space Flight Dynam- 
ics , Darmstadt, Germany, pp. 419-425, European Space Agency, 
December 1991. SP-326. 

6 M. Bello-Mora and J. Rodriguez-Canabal, “On the 3-d 
Configurations of the Cluster Mission,” ESA Symposium on 
Space Flight Dynamics , Darmstadt, Germany, pp. 471-479, Eu- 
ropean Space Agency, December 1991. S P-326. 

7 D. Hocken and J. Schoenmaekers, “Optimization of Cluster 
Constellation Manoeuvres,” 16th International Symposium on 
Space Flight Dynamics , Pasadena, California, December 2001. 

8 S. Curtis, “Auroral Lites, A MIDEX International Multi- 
probe Mission Proposal,” NASA Goddard Space Flight Center, 
Greenbelt, Maryland, August 1998. 

9 M. Hametz, D. Conway, and K. Richon, “Design of a For- 
mation of Earth Orbiting Satellites: The Auroral Lites Mission, 
Proceedings of the 1999 NASA GSFC Flight Mechanics and Es- 
timation Conference , Greenbelt, Maryland, pp. 295-308, 1999. 

10 S. Curtis, “The Magnetospheric Multiscale Mis- 
sion... Resolving Fundamental Processes in Space Plasmas,” 
NASA Goddard Space Flight Center, Greenbelt, Maryland, 
December 1999. NASA/TM — 2000-209883. 

11 A. Edery and C. Schiff, “The Double Lunar Swingby of the 
MMS Mission,” 16th International Symposium on Space Flight 
Dynamics , Pasadena, California, December 2001. 

12 T. Yang, Finite Element Structural Analysis. New Jersey: 
Prentice-Hall, Inc., 1986. 

13 P. Robert, A. Roux, C. Harvey, M. Dunlop, P. Daly, and 
K. Glassmeier, “Tetrahedron Geometric Factors,” Analysis 
Methods for Multi- Spacecraft Data (G. Paschmann and P. Daly, 
eds.), pp. 323-348, Noordwijk, The Netherlands: ISSI Report 
SR-001, ESA Publications Division, 1998. 

14 C. Russell, “Geophysical Coordinate Transformations,” 
Cosmic Electrodynamics, No. 2, 1971, pp. 184-196. 

15 C. Barton, “Revision of International Geomagnetic Refer- 
ence Field Released,” http://ww.agu.org/eos_elec/95242e. 
html, 1996. American Geophysical Union. 

16 M. Hametz, D. Conway, and D. Weidow, “New Tools for 
Formation Monitoring and Maintenance,” 1999 American Con- 
trol Conference , San Diego, California, June 2-4 1999. 

17 J. Verner, “Explicit Runge-Kutta Methods with Estimates 
of the Local Truncation Error,” Society for Industrial and 
Applied Mathematics (SIAM) Journal of Numerical Analysis , 
Vol. 15, August 1978, pp. 772-790. 

18 C. Schiff, D. Rohrbaugh, and J. Bristow, “Formation Fly- 
ing in Elliptical Orbits," 2000 IEEE Aerospace Conference, Big 
Sky, Montana, April 18-25 2000. Paper 317. 


19 L. Mailhe, C. Schiff, and S. Hughes, “Formation Flying 
in Highly Elliptical Orbits: Initializing the Formation,” Pro- 
ceedings of the International Symposium on Space Dynamics , 
Biarritz, France, CNES, June 26-30 2000. Paper MSOO/21. 

20 S. Hughes and L. Mailhe, “A Preliminary Formation Flying 
Orbit Dynamics Analysis for Leonardo-BRDF,” 2001 IEEE 
Aerospace Conference , Big Sky, Montana, March 10-17 2001. 
Paper 905. 

21 C. Chao, J. Pollard, and S. Janson, “Dynamics and Control 
of Cluster Orbits for Distributed Space Missions,” AAS/AIAA 
Spaceflight Mechanics Meeting , Breckenridge, Colorado, Febru- 
ary 7-10 1999. AAS 99-126. 

22 K. Alfriend and H. Schaub, “Dynamic and Control of 
Spacecraft Formations: Challenges and Some Solutions,” Jour- 
nal of the Astronautical Sciences, Vol. 48, April-September 
2000, pp. 249-267. 


6 


American Institute of Aeronautics and Astronautics 



FORMATION 



Fig. 3 Initializing at BLJ2 Mean Anomaly of 0 degrees 
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Fig. 7 GSM Initial Formation; TA= 180 deg 
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Formation Evolution; TA= 270, 360, 90, and back to 180 deg (left to right) 
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